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Abstract. We study the dynamics of the limit order book of liquid stocks after experiencing large intra-day 
price changes. In the data we find large variations in several microscopical measures, e.g., the volatility 
the bid-ask spread, the bid-ask imbalance, the number of queuing limit orders, the activity (number and 
volume) of limit orders placed and canceled, etc. The relaxation of the quantities is generally very slow that 
can be described by a power law of exponent ~ 0.4. We introduce a numerical model in order to understand 
the empirical results better. We find that with a zero intelligence deposition model of the order fiow the 
empirical results can be reproduced qualitatively. This suggests that the slow relaxations might not be 
results of agents' strategic behaviour. Studying the difi'erence between the exponents found empirically 
and numerically helps us to better identify the role of strategic behaviour in the phenomena. 

PACS. PACS-key discribing text of that key - PACS-key discribing text of that key 

1 Introduction a large earthquake: The number of aftershocks is described 

by a power law and no typical time scale of the relax- 



Understanding the relaxation of a system to its typical 
state after an extreme event may give a lot of information 
on the dynamics of the system. Perhaps the first study of 
power law relaxations after extreme events was by Omori 
who was studying earthquake dynamics [T]. The Omori 
law describes the temporal decay of aftershock rates after 



ation can be found. Several further examples can be found 
for non exponential relaxations in complex systems: con- 
densed state systems [5], spin glasses |3j, microfracturing 
phenomena [4 , internet traffic |5|6j . etc. In case of finan- 
cial markets Ref. showed that volatility after market 
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crashes can be also characterised by a power law decay, 
which is often referred to as the financial Omori law. 

Our aim is to analyse large price changes that can hap- 
pen often, maybe every month in case of liquid stocks. It 
is important to stress that we are not interested in market 
crashes or bubbles (when such changes happen throughout 
the market) but in large intra-day price changes specific to 
a particular stock. Understanding how the market relaxes 
after extreme events may be very important for volatility 
modeling and forecasting. 

Refs. |8l9j analysed the post-event dynamics of large 
price changes appearing on short time scales for the New 
York Stock Exchange and Nasdaq. They found sharp peaks 
in the bid-ask spread, volatility and traded volume at the 
moment of the events with slow decay to normal value 
that in some cases could be characterised by a power law. 
Several studies dealt with the analysis of the structure 
of the limit order book preceding a large price change 
|10I11I12I13| . Their results show that the volume of mar- 
ket orders play a minor role in the creation of large price 
jumps. Instead it is the disappearance of liquidity in the 
limit order book that results in extreme price changes. Ref. 
[12] also studied the relaxation of the bid-ask spread after 
large bid-ask spread variations. They found a slow relax- 
ation to normal values, characterised by a power law with 
exponent around 0.4-0.5. Ref. [T3j cross-correlate high- 
frequency time series of stock returns with different news 
feeds, showing evidence that news cannot explain the price 
jumps: In general jumps are followed by increased volatil- 
ity, while news are followed by lower volatility levels. 
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In this paper we study the dynamics of the limit order 
book of the London Stock Exchange (LSE) around large 
price changes. Analysing the limit order book allows us 
to look at the market at the level of single orders, and 
this way connect the microscopic dynamics to macroscopic 
measurables and possibly human behaviour. We focus on 
the post-event dynamics and the relaxation of the different 
measures. 

We study the dynamics of the volatility and the bid-ask 
spread near large events. We find that both have peaks at 
the moment of the price change. Their relaxation is slow 
and can be characterized by a power law. Analysing the 
behaviour of market participants, we show results on the 
bid-ask imbalance, the number of queuing limit orders in 
the book, the aggregated number and aggregated volume 
of limit orders arriving, the aggregated number of cance- 
lations and the relative rate of different order types. We 
find that the shape of the book and the relative imbalance 
changes very strongly, with a peak at the event and slow 
decay afterwards. The activity of both limit orders arriv- 
ing and being canceled increases and, after a peak at the 
event, decays according to a power law. Surprisingly we 
find that the relative rates of limit orders, market orders 
and cancelations do not vary strongly in the vicinity of 
price jumps. 

For the relaxation of most of the above measures we 
find power laws with similar exponents around 0.3 — 0.4. 
The exponents are very similar, suggesting that there might 
be a common cause behind the slow relaxations of the dif- 
ferent measures. 
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To further study the possible reasons for the relaxation 
of the volatility and the bid-ask spread, we construct a zero 
intelligence multi-agent model mimicking the actual trad- 
ing mechanism and order flow. When introducing large 
price jumps in the model, we find slow relaxations in both 
of the values. This suggests that the slow relaxations can 
be reproduced without complicated behavioural assump- 
tions. We show analytic results on relaxation of the bid-ask 
spread in the zero intelligence model. 

The paper is built up as follows: In Section [5] we re- 
view the continuous double auction mechanism and some 
properties of the limit order book. In Section [3] we present 
the data set and explain the method of determining large 
events. Section [5] shows our empirical results. In Section 
[5] we present our model and review the numerical and an- 
alytical results. We close the paper by summarising our 
results and present the conclusions in Section [Sj 

2 The limit order book 

The market we are studying is governed by a continuous 
double auction. Continuous, since orders to buy or sell can 
be introduced any time to the market and are matched by 
an electronic system in the moment when a match be- 
comes possible. In the market agents can place several 
different types of orders. These can be grouped into two 
main categories: limit orders and market orders. 

Patient traders may submit limit orders to buy or sell 
a certain amount of shares of a given stock at a price 
not worse than a given limit price. Limit orders are not 
necessarily executed at the moment they are submitted. 
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In this case they are stored in the queue of orders, the 
limit order book. 

On the other hand impatient traders may put market 
orders, orders to buy or sell a certain amount of shares 
of a given stock at the best price available. Market orders 
usually are followed by an immediate transaction, matched 
to existing limit orders on the opposite side of the book 
according to the price and the arrival time. 

The third important constituent of market dynamics 
are cancelations |14| : that is removing an existing limit 
order from the book. In general, limit orders increase liq- 
uidity, while market orders and cancelations decrease liq- 
uidity. 

It is common to analyse effective limit and market or- 
ders, i.e. regard all orders that result in an immediate 
execution as market orders. This is right from the point 
of view of the effect of these orders. However, since we 
are interested in the behaviour of traders and in the way 
their order putting strategies may change, we do not use 
this notation, and define limit and market orders by their 
code in the limit order book, thus by the intention of the 
traders. 

2.1 Notations 

Buy limit orders are generally called bids and sell limit 
orders are generally called asks. At any time instant there 
exists a buy order with highest price (highest bid), bt, and 
a sell order with lowest price (lowest ask), a,, in the limit 
order book. The mean of the best bid and ask prices is the 
mid-price: 
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We define volatility as the absolute change of the log- 
arithm of the mid-price: 

X, = |logm, -logmr_i|, (2) 

where time is measured in minutes. We apply the abso- 
lute value of the return for the non-averaged volatility in 
contrast to its square (which leads after averaging to the 
square of the standard deviation) . The difference between 
the logarithms of the lowest ask and the highest bid is 
called the bid-ask spread: 

St=\oga,-\ogb,. (3) 

This gives a measure of transaction costs in the market 
(and on the other hand the profit of market making strate- 
gies). 

Further important factors in the limit order book are 
the gaps, i.e. the number of adjacent unoccupied price 
levels between existing limit orders. Most often one talks 
about the first gap, defined as the difference between the 
best log-price and the next best log-price in either side of 
the book: 

g(" = |log;,f-'-log;,ri- (4) 

The gap is defined in the same way for both sides of the 
book and their statistics on the two sides are very similar. 
The gap is one tick if adjacent price levels are filled. The 
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second (g'^'), third (g'^'), etc. gaps can be defined in a 
straightforward way. 

3 Data and methodology 

3.1 The data set 

We studied the data of 12 liquid stocks of the LSE for 
the period 05.2000 to 12.2002. The stocks studied were: 
Astrazeneca (AZN), Baa (BAA), Boots Group (BOOT), 
British Sky Broadcasting Group (BSY), Hihon Group (HG.), 
Kelly Group Ltd. (KEL), Lloyds Tsb Group (LLOY), Pru- 
dential (PRU), Pearson (PSON), Rio Tinto (RIO), SheU 
Transport & Trading Co. (SHEL), Vodafone Group (VOD). 

The LSE consists of two parts: the downstairs market 
(SETS) and the upstairs market (SEAQ). The downstairs 
market is governed by a fully automatic order matching 
system, while the upstairs market is arranged informally 
through direct connections between agents. We confine our 
study to the electronic downstairs market. On the LSE 
there are no official market makers, instead every trader 
can act as a market maker anytime by posting bid and ask 
orders simultaneously. For more details on the rules of the 
LSE see [H]. 

3.2 What are large events? 

To study the limit order book dynamics around large events, 
we first have to define a filter that determines not only 
large price changes but also the moment of the event in a 
consistent manner. When determining the events, one has 
to face the following problem: there are volatile stocks for 
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which even a price change of 3-4 % can be an everyday 
event, while for some less volatile stocks a much smaller 
price change can be the sign of a major event. In filter- 
ing for large events, we are going to follow the method 
proposed by [9]. 

In order to determine the events, we use two filters on 
our data. 

1. Absolute filter The first filter searches for intra-day 
price changes larger than 2 % of the current price in 
time windows not longer than 120 minutes. 

2. Relative filter We measure the average intra-day volatil- 
ity pattern for the stock in the period prior to the 
event. The filter searches for intra-day price changes 
in time windows not longer than 120 minutes, exceed- 
ing 6 times the normal volatility during that period of 
the day. 

Normal volatility is defined as the average volatility for 
the same period of the day computed over the 60 preced- 
ing days. A price change is considered an event if it passes 
both of the above filters. When looking for large events 
we use transaction prices (both for returns and volatility) 
and the price change is understood to be change in the 
log-price. Furthermore we omit the first 5 minutes of the 
trading day to exclude opening effects from our measure- 
ments. We also omit the last 60 minutes of the trading 
day in order to be able to study the intra-day after event 
dynamics. The method is quite robust when altering the 
threshold values. 

To be able to localize the exact moment of the events 
we look at the earliest and shortest of the time windows 
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in which both of the thresholds have been exceeded. This 
means that, when looking for 120 minute events, if the 
price change already passes the filters in, say, 42 minutes, 
then we assume the price change has taken place in 42 
minutes and the end of the time window is set to the ear- 
liest moment by which the thresholds have been exceeded. 
The end of the time window in which the event took place 
is considered the end of the price change and from this 
point we start studying the post event dynamics. This 
way the minute of the event is exactly the end of the 
time window. 

With the filters defined above we were able to deter- 
mine 289 events for the 12 stocks in the roughly 2.5-year 
period. We found a total of 169 downward events and 120 
upward events. 

4 Empirical results 

We studied the dynamics of several measures of the or- 
der book. In Figure [1] we show the changes in the volatil- 
ity, the bid-ask spread and the number of limit orders 
placed and canceled. All plots show the dynamics for the 
60-minute pre-event and 120-minute post-event periods. 
For all events we defined time as the end of the short- 
est and earliest time windows in which the filters defined 
in Section 13.21 have been passed. Then, in case of each 
event we compared the dynamics in the period from -60 
minutes to 120 minutes to the average volatility dynam- 
ics for the same period of the day computed over the 60 
days preceding the event. In the next step we aggregated 
the dynamics over the events. This way all the plots show 
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dynamics compared to the average dynamics, daily peri- the relaxation of the same measures after the event on a 
odicities excluded. log-log scale. In order to quantify the relaxations we ap- 



■n- h / n i , ,11 • r ii 1 i-i-i 1 ply power law fits to the curves and compare the power 
h igure 1 (a) shows the dynamics of the volatility and ^ t- t- 



TT h /I \ i 1 ii 1 • r ii 1 • 1 1 IT law exponents. (Note that to study the relaxations to the 

_bigurc l(bj shows the dynamics of the bid-ask spread. In ^ ^ 



. ■ 1 J J 1 ii 1 J J J • normal value, we plot the excess variables, thus the dif- 

this case we included both downward and upward price 

ference between the actual value and the value in normal 



jumps in our sample. Figure 1(c) and 1(d) show how the 



1 f-i- 1 1 J u u J 4-u u-j -J periods.) In Table [T] we summarise the exponents of the 

number of limit orders placed has changed on the bid side / i—i t- 

and the ask side of the book in case of downward price relaxations. 



jumps. Figure 1(e) and 1(f) show how the number of limit As we can see the volatility increases very high, to 



orders canceled has changed on the bid side and the ask roughly 12 times its normal value and its decay can be 
side of the book in case of downward price jumps. The well fit by a power law with an exponent of approximately 
number of orders placed and canceled may be regarded as 0.38. The result on the power law relaxation with expo- 
measures of the activity of market participants. (In case nent 0.38 of the excess volatility after events can be com- 
of the activities, the plots for upward price jumps are very pared to results on the relaxation after crashes. Ref. [7] 
similar to those with downward jumps, so we only show shows that volatility after stock market crashes decays 
the negative events.) approximately with an exponent of 0.3, showing that the 

J2 1 ii i ii 1 ■ • -1 A n post-crash dynamics are similar to those of earthquakes, 

Wc hnd that the dynamics are very similar. All mea- ^ j j 

1 , liii iTi-i commonly known as the Omori law. Another study showed 
sures have a strong peak at the moment oi the event deter- 

J u cii- • i r ii 1 i- that the relaxation of excess volatility after fluctuations is 
mined by our hitcr. We are going to focus on the relaxation 

r.i . 1 1 n ;i 1 rpi ■ • ,1 i-i. ■ n characterized by a power law with exponent between 0. 3- 

of the variables after the peak: i his IS the part that IS well j f t- 

„ 11 izii • ii 1 /-I • ii ■ • 0.4, measured to be robust for different time periods and 
dehned by our hitermg method. Concerning the rise prior 

. ,1 , 1 ,1 r 1 1 T-v i across markets 1161 . This may point us to the fact that 

to the events we have to be careful, however. Due to our ' — ' ^ 

^11^1, . . i r i ii i 1 i2 ii relaxation of the volatility after extreme events and after 
method of determining events (the fact that we dehne the 

1 ^,1 ^. .1 ^. r ^^ i^ ii fluctuations are similar, also mentioned in 1171 . 

end of the time window as zero time of the event) the pre- ' — ' 

event dynamics are conditioned on the event and can not The bid-ask spread increases, but to a lower value than 

be regarded as independent. Practically, what we can say the volatility, with a peak of roughly 3 times the average. 

is that the variables change near the event with a peak at The first part of the relaxation curve seems to show scaling 

the moment of the event, but the actual rise can not be behaviour and may be fit by a power law. Interestingly, we 

characterised through these results. The relaxations after find that the exponent of the decay is around 0.38, very 

the event seem to be slow in all cases. In Figure[5]we show close to the exponent of the volatility decay. The slow. 
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Fig. 1. (Color online) Dynamics around events compared to the usual dynamics. Usual dynamics are computed from the 60-day 
pre-event average for the same period of the day. | (a) | Mid-price volatility averaged over all events; [(b)] Bid-ask spread averaged 
over all events; [(c)] The number of limit orders in case of negative events on the bid side of the book and |(d)| on the ask side of 
the book; |(e)| The number of limit order cancelations in case of negative events on the bid side of the book and |(f)| on the ask 
side of the book. 
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Fig. 2. (Color online) Relaxation of the excess variables on log-log scale with power law fits: |(a)| Mid-price volatility averaged 
over all events, exponent; 0.38 ±0.01; |(b)| Bid-ask spread averaged over all events. The first part of the relaxation curve seems 
to show scaling behaviour and may be fit by a power law with an exponent: 0.38 ±0.03. However the scaling for the bid-ask 
spread is less consistent than that of the volatility. Also the long time decay is apparently even slower than a power law.; WE 



The number of limit orders in case of negative events on the bid side of the book, exponent; 0.37 ±0.01, and (d) on the ask side 
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power law like decay of excess bid-ask spread is in agree- 
ment with the findings of [T^ . However the scaling for the 
bid-ask spread is less consistent than that of the volatility. 
Also the long time decay is apparently even slower than 
a power law. Interestingly, [9] found strong variation of 
the bid-ask spread after large events at NYSE, but not 
on Nasdaq. They argued that this difference is caused by 
the diverse trading rules on the two markets, pointing out 
that the existence of a single market maker at NYSE leads 
to the opening of the bid-ask spread. The trading rules of 
LSE, being a fully automatic market makes it similar to 
Nasdaq, so, according to that argument, one would ex- 
pect low variation of the bid-ask spread on the LSE. This 
is not consistent with our findings. An explanation of this 
contradicting behaviour on the two markets is missing at 
present. 

The number of limit orders placed increases both on 
the bid and the ask side of the book. We find a peak 
of roughly four times the usual value with a slow decay 
of the activity afterwards. The relaxation of the excess 
limit order putting activity after events can be described 
by a power law decay with an exponent of roughly 0.37- 
0.4. Concerning the dynamics of the aggregated volume of 
limit orders placed we get very similar results (not shown 
here). However, in case of the volumes we see a slightly 
stronger increase around the event with roughly 6 times 
the usual value. This suggests that the average volume 
of limit orders increases around large price changes. The 
relaxation of the excess limit order volume can be charac- 



order book around large price changes 9 

Table 1. The exponents of the relaxation curves for the dif- 
ferent variables. 



variable 




exponent 


volatility 




0.38±0.01 


bid-ask spread 




0.38±0.03 


limit orders placed - 


bid 


0.37±0.01 


limit orders placed - 


ask 


0.40 ±0.01 


cancelations - bid 




0.30 ±0.02 


cancelations - ask 




0.42 ±0.02 



terised by a power law of exponent roughly 0.44 in case of 
negative events and exponent 0.48 for positive events. 

Similarly, the number of limit orders canceled increases 
strongly around events, with a slow decay after the peak. 
It also seems that the increase is stronger and much more 
clear on the side of the book which is opposite to the di- 
rection of the price change, i.e. ask side in case of negative 
events and bid side in case of positive events. For the relax- 
ations we find that in general, the decays on the opposite 
side of the book compared to the direction of the price 
change can be well fit by a power law of exponent roughly 
0.42. The decays on the same side as the price change di- 
rection are more noisy but may be fit by a power law, 
showing exponents of roughly 0.3-0.35. 
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It is generally understood that the relative amount of show that there is a huge imbalance of volume around 

supply and demand govern the movement of prices. To events. The relaxation of the imbalance after the events is 

quantify the pressure of orders from either side of the very slow. 

book, we study the bid-ask imbalance on the market. We ^ measure very similar to the bid-ask imbalance, is the 

denote the total volume to buy on the market by y/'"' and riumher of queuing orders on either side of the limit order 

/sell 



the total volume to sell by Vf' . Then we define the buy ^ook. Figure lOT shows the dynamics of the number of 



imbalance and the sell imbalance as: queuing orders in the limit order book compared to the 

60 day pre-event interval for the same period of the day, 

jbuy /t-\ 

' -i/buy .^rseii ^ ' for negative events. The plot shows results for both the 

and bid and the ask sides of the book. The results for positive 

events are very similar (of course symmetrically to the 

-ysell 

f^" — t (Q) case of negative events), so we only show the figures for 

' y*">' -I- ' 

Trivially, ^ ^ the negative events. 

In Figure [3] we show the dynamics of the imbalances We can see that for negative price jumps the number 

for upward and downward price changes, compared to the of queuing limit orders on the bid side decreases to about 

60 day pre-event interval for the same period of the day, ^^^^^ usual value and only relaxes back slowly. At the 

as before. The two curves are very similar showing a van- same time on the ask side the number of queuing orders 

ishing amount of orders to sell in case of upward price increases to roughly 1.6 times the usual value and even 

moves and vanishing amount of orders to buy in case of after the event stays very high for a long time. Further 

downward price moves studies are needed to see if this very slow relaxation is a 

To understand the absolute values on the y-axis is not sign of the limit order book being partly frozen in post- 
straightforward, to make it clear we show the following: If event peiiods. 

we assume that in regular market periods, on average half It is interesting to study the rates of different market 
of the volume of orders is to buy and consequently half is activities, i.e. the relative number of limit orders, mar- 
to sell, then the values on the y-axis show that in case of ket orders and cancelations, compared to the total num- 
the moment of a positive event, only w 100% — 50%* 1.57 = ber of orders. Figure [4(b)| presents the dynamics of the 
21.5% of the total volume of orders appears on the sell side three rates around negative events for the buy side of the 
of the book and in case of the moment of a negative event, book. We see that there are no strong changes in the rate 
only w 100% — 50%* 1.48 — 26% of the total volume of around large price changes. The results for the ask side 
orders appears on the buy side of the book. These numbers of the book are the same. When studying positive events 
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Fig. 3. (Color online) The dynamics of the imbalance in the volume of supply and demand: [(a)] shows the buy imbalance (/*"-' ) 
in case of upward price jumps; (b) shows the sell imbalance (/'''") in case of downward price jumps. 
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Fig. 4. (Color online) |(a)| Dynamics of the number of queuing orders in the book in case of negative events. The number of 
queuing limit orders on the bid side decreases to about half of the usual value and only relaxes back slowly. The number of 
queuing limit orders on the ask side increases to roughly 1.6 times the usual value and even after the event stays very high 



for a long time.; (b) The dynamics of the rates of limit orders (black circles), market orders (red squares) and cancelations 



(blue triangles) around negative events for the buy side of the book. No very strong variations can be seen in the different rates 
around the large price changes. 



we get very similar dynamics: No strong variation in the of large events and not the strategy of traders how they 
relative rates. This result means that it is rather the entire place ordersO 
market activity changing (increasing) in the surroundings 



^ As we stated before, we regard orders as limit or market 
orders by the intention of the trader, not by their effect (i.e. 
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Summarising the empirical results: We have studied 
the dynamics of several measures of the order book be- 
fore and after large price changes. For the change in the 
volatility, the bid-ask spread, the limit order placing and 
cancelation activity, the bid-ask imbalance, the number of 
queuing orders in the book we found strong variation at 
the moment of the event and slow relaxation in the post- 
event period. Specially, in case of the volatility, the bid-ask 
spread and the activities we found a relaxation very simi- 
lar to a power law, with exponents close to 0.4 suggesting 
a possible common cause behind the slow relaxations of 
the different measures. 

Analysing the rates of limit orders, market orders and 
cancelations around large events, we did not find strong 
variation, showing that strategy of traders in choosing 
their type of order does not vary much. 

5 An agent-based model 

As stated above, we found very similar relaxation in dif- 
ferent measures of the limit order book. To better under- 
stand the dynamics leading to the slow relaxations, in this 
section we introduce a multi-agent model of the order plac- 
ing and removing process. When constructing a modeling 

framework, we have to decide which path to follow: 
not effective limit and market orders). Because of this we get 
a higher rate of limit orders and lower rate of market orders, 
than presented in 10 . However the rate of cancelations fits 
their results very well, so we believe that after accounting for 
the effective orders, there are no contradictions between the 
two results. 



order book around large price changes 

1. Building a multi- agent model with complicated strate- 
gics, involving behavioural assumptions. 

2. Building a zero intelligence multi-agent model, with 
some very basic assumptions on the order flow. 

In the literature there are several examples for both types 
of models. Models of type 1 permit one to study behavioural 
results of the model, but when building the trading strate- 
gies we have to be careful, not to assume unrealistic prop- 
erties of traders and/or avoid the common error, to simply 
find as output exactly the input assumptions. Models of 
type 2 are easier to construct, but apart from the problem 
of possible over-simplification, they also confine us to the 
analysis of non-behavioural measures through the model. 

We chose to follow the path of |l8ll9l20l2ll22l23l24l25| 
to construct a zero intelligence multi agent model of the 
continuous double auction through the limit order book. 
The model is aimed to be as simple as possible but captur- 
ing the most important properties of the continuous dou- 
ble auction. If we are able to reproduce some results with 
a zero intelligence model it may suggest that the particu- 
lar phenomenon is not due to traders' strategic behaviour 
but rather to the market mechanism or institutions. 

5.1 Details of the model 

We assume limit order placing and cancelations similarly 
to a deposition-evaporation process and furthermore in- 
troduce market orders. Our model is similar to Maslov's 
model [5D]. The main differences are that we allow for 
cancelation of existing limit orders (and through this the 
tuning of the probabilities of different actions) and that 



Bence Toth et al.: Studies of the limit order book around large price changes 



13 



evaporation = 
order cancellation 



market order: 
eliminating 
opposite best 
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deposition = limit 
order arriving 



agents put their limit orders relative to the mid-price, this mechanism. Since in our model there are no crossing limit 

way we allow for a non trivial dynamics of the bid- ask orders, in order to have a stationary number of orders we 

spread. All orders arrive or evaporate with the same unit set P^o = 0.5. 
volume. The trading mechanism is the following: 

— Limit orders arrive with rate Plo per unit time with 
equal probability to buy or sell. Limit orders get de- 
posited in the interval [m, — D,m,] in case of buy orders 
and in the interval [nit , nit +£^1 in case of sell orders with 
uniform distribution, where nit is the mid-quote price 
(see Equation [T]) and D is a parameter of the model. 

— Market orders arrive with rate Pmo per unit time with 
equal probability to 'buy' or 'sell'. A market order 
to buy (sell) will get executed immediately by being 
matched to the best limit order to sell (buy). 

— Existing limit orders are being canceled with rate Pc 
per unit time from the 'buy' or 'sell' side with equal 
probability. In case of a cancelation on one side of the 
book, all limit orders on that side have the same prob- 



m(t) 
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ability: p = V,J^,),/ to be canceled where Vtotai is the total 
volume of limit orders on that side of the book. Thus 
on average one limit order evaporates from the book 
in case of cancelations. 

The three rates add up to one: Plo + Pmo + Pc = 1 (in 
other words we study the model in event time). In the 
beginning of the simulation there is a long warming up 
period, when orders are randomly placed in the book, in 
order not to have spurious results due to fluctuations and 
an empty book. Similarly to the empirical analysis, we 
denote the mid-price at time t by m, and the bid-ask spread 
at time f by 5,. Figure[5]shows the scheme of the order flow 



Fig. 5. (Color online) Scheme of the dynamics of the zero 
intelligence model. 

We generate large price changes manually, by clearing 
out all limit orders on one side of the booj^ in the interval 
[bt—J,bt\ for price drops and [at^Ut+J] for price jumps, 
where 7 is a parameter of the model (and bt and a, are the 



^ This method is clearly not realistic. However, this is the 
only way we could think to generate large price changes in a 
zero intelligence model, i.e., when price jumps are not governed 
by sudden changes in the agents perception, strategy and herd- 
ing behaviour. Since we are interested in the relaxation after 
the jump, we believe that this way of generating the large price 
changes is satisfactory. 



14 Bence Toth et al.: Studies of the limit 

best bid and best ask respectively), this way the jump in 
the mid-price is 7/2. 

As we have mentioned, with such a simple model we 
have to confine ourselves to the study of the bid-ask spread 
and the volatility, i.e., variables whose dynamics may be 
studied in a model without adaptation rules and agents' 
intelligence (unlike other variables, like the activity, that 
are more closely related to changes in agents' perception 
and thus cannot be analysed in a simple zero intelligence 
model). 

5.2 Numerical results 

For usual market dynamics the rate of different orders are 
roughly Plo = 0.5, Pmo = 0-16 and Pc 0.34 loE In this 
section we present the results of our numerical simulations 
using the above empirical probabilities. The parameters of 
the simulation are the following: 

- D= 1000 

- / = 1000 

— the frequency of large price jumps was / = (5 * lO''^)^' 

— the length of the simulation was 5*10"^. 

Figure [S] shows the average decay of the volatility and 
the bid-ask spread after large price changes in the model. 
The plots are created in the manner as the figures in Sec- 

Note that these rates are defined for effective limit orders 
and effective market orders, i.e. all orders that lead to an im- 
mediate execution are regarded as market orders. Since we do 
not have crossing limit orders in our simulation, it is right to 
use these values. 
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tion 131 Time zero is the moment of the event and the y- 
axis shows the relative dynamics compared to stationary 
market periods, averaged over 100 events. The decays are 
qualitatively similar to what we have seen for the empir- 
ical data. The short time relaxations (up to roughly 100 
simulation steps) can be described by power laws, with 
exponents of roughly 0.5 (0.50 ±0.04 for the volatility and 
0.48 ±0.01 for the bid-ask spread). However, there is a dif- 
ference between the empirical and numerical exponents. 
This suggests that our simple model is not able to en- 
tirely reproduce the relaxations. It seems that part of the 
slow relaxation can be generated in the model, however 
the difference between the exponents 0.4 and 0.5 is im- 
portant and we believe that the discrepancy is due to the 
behaviour of the agents, that can not be captured in the 
zero intelligence model. It is interesting that the exponent 
of the decay in the volatility and in the bid-ask spread 
seem to be very close to each other, similarly as in the 
case of the empirical data. Also similarly to the empirical 
data, the peak in the volatility is smaller than that in the 
bid-ask spread. 

5.3 Analytical treatment 

As we have seen, the bid-ask spread shows a slow relax- 
ation both for empirical data and simulations. To under- 
stand if it is really a critical relaxation and what deter- 
mines the exponent, we tried to treat the model analyti- 
cally. 
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D=1000; J=1000 



D=1000; J=1000 




1000 




Fig. 6. (Color online) Relaxations after price jumps in the numerical model, in case of Pio ~ 0.5, Pmo = 0.16 and P(j — 0.34. We 



also show the power law fits of the relaxations, (a) relaxation of the volatility, 0.50 ±0.04; (b) relaxation of the bid-ask spread. 



0.48 ±0.01; Both exponents are close to 0.5. Similarly to empirical results, the variation in the volatility is much stronger than 
in the bid-ask spread. 



5.3.1 A limit case 



bid-ask spread. With the above assumption, the expecta- 
tion of the change in bid-ask spread from one time step to 
another can be written in the following way: 



k=\ 



(7) 



The model can be treated analytically in the case of Pmo = 
0, i.e. when only limit order placing and cancelation de- 
termine the flow {Plo = Pc = 0.5), using some simple as- 
sumptions. 

Mathematically the bid-ask spread in this limit case 
can be understood as the time evolution of the minimum ^^ere the first term on the right hand side assumes that 
of numbers distributed uniformly on a finite interval (ex- cancelations do not change the bid-ask spread, and the sec- 
cept for very short times). This tells us that the bid-ask ^nd term describes the effect of a new limit order arriving: 
spread will decay according to a power law with unit ex- jf falls outside the bid-ask spread it does not change the 
ponent, clearly showing that the limit case is unrealistic, bid-ask spread, if it falls inside the bid-ask spread with a 
Nevertheless we discuss it briefly as we are going to use distance k from the same best price, it decreases the bid- 



these results in the following. 

When describing the time evolution of the bid-ask spread j^] -^q gQ\^ 
we use a simplifying assumption: Since the probability of 
the best order to be canceled on either side is very low, we 
can assume that cancelations do not alter the value of the 



ask spread exactly by k. Summing up the right side of Eq. 



= -P,j'l^^\_ 



E(A5,) = -Plo 
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Using Equation \8\ in case of Plq = 0.5 the expected account for the gaps as well. The expectation for the gap 
value of the bid-ask spread can be written in the following can be given by the following equation: 



recursive formula: 



E(5,+i)=E(5,) + E(A5,) = 
^ 8dJ 16D 



(9) 



The above recursive formula is a mean field theory that 
can be used to describe the relaxation when knowing the 
initial bid- ask spread after the event. The formula fits the 
numerical results very well reproducing the asymptotic 
power law with exponent very close to 1. 

5.3.2 General case 

We are most interested in the general case when the prob- 
ability of market orders is finite. For Pmo 0, the change 
of the bid-ask spread is the following: 



k=\ 



+Pmo — ^Plo 



W AD 



PMO-gY\ (10) 



where the first term assumes again that cancelations do 
not change the bid-ask spread, the second term gives the 
expected change in case of a limit order arriving (simi- 
larly to the limit case, Eq. [7]) and the last term stands 
for a market order arriving, increasing the bid-ask spread 
exactly by the size of the first gap, (the gap is defined 
in the same way as in the empirical case). As we can see, 
when having market orders, not surprisingly we have to 



+Plo 



D 



D 



-yk+ 

D H 



k=l 



k=l 



- Pmo ■ sl^i , 



(11) 



where g^^^ stands for the second gap in the book, i.e. the 
price difference between the second best and third best or- 
der on the same side of the book. In Equations [TU] and [11] 
we neglected the probability of cancelations changing the 
bid-ask spread or the first gap. As we can see, in the ex- 
pected value of the first gap, a term containing the second 
gap occurs. This is the general case for all gaps that is, 
when writing up g^"\ it will contain a term depending on 
gi"+^)_ ■\Ye have to find a closure for this infinite hierarchy 
of equations, so we do the following. First, we estimate 
the relation between the first and the second gaps for the 
equilibrium state of the system (the value of the second 
gap, as a function of the first gap). Second, we assume this 
relation to be constant also for the relaxation period. By 
equilibrium state, here we mean the long time behaviour 
of the system without large price jumps. We denote the 
equilibrium values of the bid-ask spread, the first gap and 
the second gap by o, and "/^^ respectively. In case of 
stationarity, E(A5';) =0, thus Equation [TUl becomes: 



Plo 

Equation [TT] becomes 



a a 



Pmo-T 



(12) 
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+Plo 



D 



k=\ 



1 ^ 



PMO-i 



12) 



(13) 



Combining Equations [T^ and [T51 we get the following for- 
mula: 



Y(i) IDPmo 



2Pw 



(14) 



Knowing that Plo = 0.5 in the model, we get 



v(2) 



1 Plo 



o + Y 



(15) 



Y*') 2DPmo 

The above relation between the first and second gap 
can be assumed to be true for the relaxation process, since 
the second gap has only a minor role in case of decreasing 
bid-ask spread. Introducing Equation[TS]into Equation llll 
we get a recursive formula for the size of the first gap, 
and through that for the size of the bid-ask spread. To be 
able to use the relations, we need to know the expected 
value of the stationary bid-ask spread, o. This wc know 
from the decay of the bid-ask spread: The value, of the 
first point (the relative opening of the bid-ask spread) is 
approximately (o + 7) /o. 

Figure [7] shows the comparison of the numerical and 
the analytical results for the relaxation of the bid-ask 
spread for Pmo = 0.16. 

The analytical formula seems to describe the relax- 
ation of the bid-ask spread for short times but not long 




Fig. 7. (Color online) The relaxation of the bid-ask spread 
after large price changes in case of P^^q = 0. 16. The black circles 
show the numerical values, the red line shows the analytical 
values. 

times. We can see that the agreement between the analyt- 
ical and numerical curves only hold for roughly the first 
100 simulation steps in the relaxation. This is due to the 
fact that on long times the "mean field" assumption ig- 
noring the actual dynamics of the second gap gives rise to 
larger errors. 

6 Conclusions 

We have studied the dynamics of several measures of the 
order book after large price changes. For the change in the 
volatility, the bid-ask spread, the limit order placing and 
cancelation activity, the bid-ask imbalance, the number 
of queuing orders in the book we found strong variation 
at the moment of the event and slow relaxation in the 
post-event period. Specially, in case of the volatility, the 
bid-ask spread and the activities we found a relaxation 
very similar to a power law, with exponents close to 0.4. 

Analysing the rates of limit orders, market orders and 
cancelations around large events, we did not find strong 
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